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We examine the changeover in the particle configurations and the dynamics in dense 
Lennard-Jones binary mixtures composed of smaU and large particles. By varying the com- 
position at a low temperature, we realize crystal with defects, polycrystal with small grains, 
and glass with various degrees of disorder. In particular, we show configurations where small 
crystalline regions composed of the majority species are enclosed by percolated amorphous 
layers composed of the two species. We visualize the dynamics of configuration changes using 
the method of bond breakage and following the particle displacements. In quiescent jammed 
states, the dynamics is severely slowed down and is highly heterogeneous at any composi- 
tions. In shear, plastic deformations multiply occur in relatively fragile regions, growing into 
large-scale shear bands where the strain is highly localized. Such bands appear on short time 
scales and change on long time scales with finite life times. 

§1. Introduction 

Binary particle systems with size dispersity exhibit complicated phase behavior 
depending on the temperature T, the average number density n = {Ni -\- N2)/V , 
and the composition c = N2/{Ni -\- N2)i^^ Here A^i and are the numbers 
of the small and large particles in a volume V , respectively. At high densities, 
the particle configurations sensitively depend on the size ratio 02! (y\ between the 
diameters of the two species, o\ and a^- If (^2/(^1 is close to unity at large n and at 
low T, a crystal state is realized with a small number of defects. However, if ct2/o"i 
considerably deviates from unity, crystal states are realized only for very small c or 
1— c. For not small c and 1— c, polycrystal and glass states emerge without long-range 
crystalline order at low T. It is then of great interest how the particle configurations 
and the dynamics change with varying the composition c at considerably large size 
dispersity. As the glass transition is approached, the structural relaxation time 
Ta{T) grows dramatically from a microscopic to macroscopic time, while the particle 
configurations remain random yielding the structure factors similar to those in liquid. 

The dynamics of supercooled liquid and glass have been studied extensively 
using molecular dynamics simulations (MD). As a marked feature, the glass dynam- 
ics is highly heterogeneous EUl^ In particular, Yamamoto and one of the present 
authorsHf examined breakage of appropriately defined bonds. The broken bonds 
accumulated in long time intervals are analogous to the critical fiuctuations in Ising 
systems such that their structure factor may be fitted to the Ornstein-Zernike form, 

5b(fc) = 5b(0)/(l + fcV), (1) 

in two dimensions (2D) and in three dimensions (3D). The wave number k is smaller 
than the inverse particle size. The correlation length ^ can thus be determined. 
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which grows with lowering T. Bond breakage events tend to take place repeatedly in 
relatively fragile regions, leading to aggregation of broken bonds on long time scales. 
Kob et a/.^^' pointed out relevance of stringlike clusters of mobile particles whose 
lengths increase at low T. In addition, the diffusion constant in glassy materials 
is strongly affected by the dynamic heterogeneitylHHIll' Afterwards, some authors 
have claimed the presence of correlation between the structural heterogeneity in 
the particle configurations and the dynamic heterogeneity on long time scalesiHIil^ 
Recently significant heterogeneity has been found in the elastic moduli in glassjlSJ'Hl) 
which is the origin of nonaffine elastic displacements for very small strains. 

Glass dynamics under shear flow has also been studied by many authors )lll'~tl3 
where the shear rate 7 is a new parameter representing the degree of nonequilibrium. 
Similar jamming rheology has been found in foam and microemulsion systems, colloid 
suspensions, and granular materialsl^^'""^^' In jammed states of molecular glass, each 
particle undergoes shear- induced configuration changes on the time scale of_7~^. At 
a low temperature T, the total rate of the configuration change is given bjHH''!!!!' 

r„(r,7)-i =r„(r)-i + Aa7. (2) 

In the right hand side, the first term is the thermal activation rate and the second 
one is the shear-induced rate, where the coefficient is of order unity. The steady- 
state viscosity r/(r, 7) is proportional to rQ,(T, 7). The following djmamic scaling 
relation for the correlation length was in accord with the simulationjl^J'IHI 

e(T,7)~T,(T,7)V^ (3) 

including the sheared case. The dynamic exponent z was estimated as 4 in 2D and 
2 in 3D. In glass a Weissenberg number may be defined by 

Wi = T„(r)7. (4) 

A nonlinear response regime emerges for Wi> 1. Around the glass transition, Tq,(T) 
grows and the regime WiS> 1 is realized even for extremely small shear, where 
the configuration change is mostly induced by shear and t] ~ 7~^. Indeed in this 
condition, shear-thinning behavior was measured in glass under uniaxial stress JlZl'll^ 
Also for near-critical and complex fluids, the criterion of nonlinear shear effects is 
given by Wi= T7 > 1 with an appropriate long relaxation time rMf Moreover, a 
number of MD simulations of sheared glass have realized self-organization of "shear 
bands" with high strain localization extending througout the system. In sinrple shear 
flow such bands are nearly along the flow or the velocity-gradient direction 
while under uniaxial stress they make an angle of 7r/4 with respect to the uniaxial 
axis.P^^^ In metallurgy, shear bands have been observed in amorphous solids under 
uniaxial stretching or compression above a yield stressil^l^llL' 

Plastic deformations in crystal and polycrystal are extremely complex and are 
still poorly understood despite their extensive research lllUllf In crystal with low- 
density defects, dislocation motions play a major role in plasticity, where they yield 
slip planes with various sizes as observed in acoustic emission experimental^f and 
by transmission electron microscopy HHI In polycrystal, the particles in the vicinity 
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of grain boundaries are relatively mobile compared to those within the grains, so 
their collective motions give rise to sliding of the grain boundaries under applied 
stressHfH^^f In these systems, plastic events take place as bursts or avalanches 
spanning wide ranges of space and time scales. In this manner, they intermittently 
release the elastic energy stored in the crystalline regions at high strain. On the 
other hand, in glass under stress, plastic events have been assumed to be spatially 
localized due to the structural disorder in glass 113^111' but our recent 2D simulation 
has shown that they frequently take place over wide areas in short timesHll In 
glass, slip elements (slip lines in 2D) do not much exceed the particle size, but they 
successively appear in their neigborhood forming large-scale aggregates. When glass 
rheology is studied in simulation, such avalanches manifest themselves in large stress 
drops in the stress-strain curve. The irregularity of the average stress versus time has 
been conspicuous in numerous simulations in the literature, even though the curve 
should become smoother with increasing the system size. Molecular glassy materials 
behave as elastic bodies on short time scales with mesoscopically inhomogeneous 
elastic modulil22l''[13 Thus large-scale release of the elastic energy at high stress is a 
common feature of plasticity in crystal, polycystal, and glass. 

In this paper, we will investigate the heterogeneity dynamics in high-density 
and low-temperature 2D model binary mixtures without and with applied shear for 
various compositions on the basis of our previous workj^'llll'll^} We are interested 
in the relationship between the dynamic heterogeneity on long time scales and the 
structural heterogeneity in the particle configurations. In a separa te p aper we have 
examined rheology and plastic deformations on various time scales'^^' 

The organization of this paper is as follows. In Sec. II, our model and our simula- 
tion method will be explained. In Sec. Ill, numerical results will be presented on the 
particle configurations and the heterogeneous dynamics for various c without shear. 
In Sec. IV, we will examine the collective dynamics under shear, which are even more 
heterogeneous than in quiescent states. 

§2. Background of Simulation 

2.1. Model 

Our two-dimensional (2D) binary mixture^ 'llll'llll' consist of two species 1 and 
2 interacting via truncated Lennard-Jones (LJ) potentials, 

- C«/3 (5) 

which are characterized by the energy e and the interaction lengths = (cq,+(T^) /2 
(a, /3 = 1,2). The r = \rj — r^l is the particle distance. Hereafter rj denote 
the particle positions. The diameters of the two species are ai and a2 and their 
ratio is fixed at a^jox = 1.4. For r > rcut = 3.2cri, we set Vai3 = and the 
constant Ca/B ensures the continuity of Vai3 at r = rcut- The particle number is 
A'^ = A'^i + = 9000 and the composition c = is varied. The system is in the 

region —L/2 < x,y < L/2. The volume is chosen such that the volume fraction of 
the soft-core regions is fixed as {Niaf + N2(T2)/L'^ = 1. We integrated the equations 
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of motion using the leapfrog method, where the time step is 0.002r with 

T = ai^/mi/e. (6) 

The mass ratio is fixed at mi/m2 = (o"i/(72)^. The following data were taken at 
T = 0.2e/kB- The space coordinates x and y, the time t, the shear rate 7, and the 
temperature T will be measured in units of cri, r, r~^, and e/kB, respectively. 

There was no tendency of phase separation in our simulation. For the parameters 
we adopt in this paper, the structural relaxation time Ta{T,c) without shear is of 
order 10'^ in glass, is longer in polycrystal, and tends to infinity in crystal from 
the decay of the self-time-correlation functionllfl The grain boundary motions are 
severely slowed down in the presence of large size dispersity, while they move rather 
fast in one-component systems. Thus Tq is longer in polycrystal than in glass. 

2.2. Orientation angle and disorder variable 

In the particle configurations at high density, a large fraction of the particles are 
enclosed by six particles. The local crystalline order may then be represented by a 
sixfold orientation® We define an orientation angle oj in the range — 7r/6 < aj < 
7r/6 for each particle j £ a using the complex number jU'lH!''!!!!' 

J2 exp[6i0,fc] = |.P',|e6^°^ (7) 

fc 6 bonded 

where the summation is over the "bonded" particle k £ /3 satisfying r = | — | < 
1.5(Tq,/3. The 9jk is the angle of the relative vector Vj — Vk with respect to the x axis. 

We also introduce another non-negative-definite variable, called the disorder 
variable, which represents the degree of the deviation from the hexagonal order for 
each particles j by 

Dj= |e6i". -e^^^'^p = 2 ^ [1 - cos 6(0^ - a^)]. (8) 

fcSbonded fcgbondcd 




c 



Fig. 1. Average disorder variables _D vs c 
without shear and with shear 7 = 10"*, 
where T = 0.2 and = 9000. 



Here Dj are in the range 15-20 for par- 
ticles around defects, while Dj = for a 
perfect crystal at low T. The degree of 
the overall disorder may be represented 
by the average over the particlesjH' 

N 

D = Y,Dj/N, (9) 

which is large in glass and liquid and 
is small in crystal. In Fig. 1, we show 
the average D versus c at T = 0.2 with- 
out shear and with shear 7 = 10~^, 
where the two curves nearly coincide. 
For this shear rate (much smaller than 
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the molecular frequency r^^), there is no proliferation of defects induced by shear, 
while the defect structure is continuously deformed by shear. The maximum of D is 
attained at c ~ 0.5. See Figs. 8 and 9 of Ref. 6 for Z) as a function of T and o"2/cji. 

The angle variable aj was originally introduced to study the thermal fluctuations 
of the hexagonal lattice structure around the melting transition in 2D& We have 
recently studied 2D melting numerically with the aid of visualization of aj and Dj j£3 
where polycrystalline patterns are apparent at the melting not in accord with the 
original theorj@' (see the end of the last section for more discussions). 

2.3. Bond breakage 

The jamming dynamics can be conveniently visualized if use is made of the bond 
breakage For each particle configuration at a time t, a pair of particles i G a 
and j G /5 is considered to be bonded if 

rij{t) = \ri{t) - rj{t)\ < A^a^p, (10) 

where aa/B = (tq, + cr^)/2. We set Ai = 1.2; then, AiUap is slightly larger than the 
peak distance of the pair correlation functions Qapif)- After a time interval At, the 
bond is regarded to be broken if 

rij{t + At)>A2aap, (11) 

where we set A2 = 1.5. In the following figures, bonds broken during a time interval 
[^1,^2] will be marked by x at the middle point Rij = ^{ri{t2) + rj{t2)) of the two 
particle positions at the terminal time t2- The structure factor of these middle points 
is written as Sb{k) in Eq.(l). 

In the original simulation using the soft-core potentialP'Hl the bond breakage 
rate was rfc(r, 7)"^ = 0.lTa{T,'y)~^ , where Ta(T,'y) is the decay time of the self- 
correlation function and behaves as in Eq.(2) in shear. The r;,(T, 7) represents the 
average bond life time and the broken bond number ANi, of the whole system is of 
order NTi,{T,'y)~^At in time interval At. For Wi> 1, the bond breakage is mostly 
induced by shear and the time average of the broken bond number per particle 
ANh/N is of the order of the average strain 'jAt in the plastic flow regime. Around 
this mean value, AN^/N increases significantly on occurrences of large-scale plastic 
events and decreases in "elastic periods" without themll3 

2.4. Averaged velocity 

Shear bands have been observed widely in complex fluids, glassy fluids, and gran- 
ular materials. To examine such strain localization in shear flow, we may introduce 
an averaged velocity Vx{y, t) as follows. We first integrate the x component of the 
momentum density Jx{x,y,t) in the flow direction: 

1 f^/^ 1 
Jx{y,t) = Y dxJxix,y,t) = -y2mjXj{t)5{y - yj{t)). (12) 

^ J -L/2 ^ ^ 

We next smooth Jxiy,t) over space and time intervals with widths Ay and At as 

-, ry+Ay/2 rt+At 
Vx{y,t) = ——^l dy' dt'Uy',t'), (13) 

pAyAt Jy-Ay/2 Jt 
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where p is the average mass density. We fix Ay at L/20. However, we choose At 
depending on the time scale of plastic events under consideration. This is because 
of the hierarchical dynamics of plastic deformations. 

§3. Numerical results in quiescent states 



Without applied shear, we first integrated the equations of motion under the 
periodic boundary conditions in the x and y directions. We equilibrated the system 
at T = 0.2 for a time interval of 5 x 10'^ until no appreciable time evolution was 
detected in various thermodynami c q uantities. At this time we will set t = 0. We 
attached a Nose-Hoover thermostal©'!^ to all the particles, which then obeyed 



_d_ 



-U — Crrijrj, 



(14) 



where fj = drj/dt, rj = d'^Vj/dt^, and U is the total potential. The thermostat 
variable C,{t) obeyed 

^ 1 r 1 ^ . 1 

(15) 



dt 



c 
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where tnh is the thermostat characteristic time. We set tnh = 0.961. 

By varying c in the range c < 0.5 at (J^jox = 1.4, Hamanaka and one of 
the present author^ 'HI realized crystal states with small numbers of defects for 




Fig. 2. (color online) Left: disorder variables Dj{t) at t = lO** with the gradation bar at the top, 
where 7 = and c — 0.9. Crystalline regions composed of the large particles (white) are enclosed 
by percolated amorphous layers containing the small particles (gray or black). Right: broken 
bonds (x) and displacement vectors for mobile particles with \Arj\ > 0.2 in the time interval 
[10*, 4 X 10*] with width At = 3 x 10* in the same run. The particles with Dj > 4 are written 
in black (left). The whole system — L/2 < x,y < L/2 is shown. 
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Fig. 3. (color online) Upper plates: disorder variables Dj{t) aX t — 10'' for c = 0.05, 0.1, and 0.15 
with 7 = 0. Lower plates: broken bonds (x) and displacement vectors for mobile particles with 
\Arj I > 0.2 in the time interval [10*, 4 x 10*] in the same run producing the corresponding upper 
panel. These are snapshots of 1/16 of the total system. The numbers at the corners (upper 
plates) denote the normalized coordinate x/L or y/L. 



c < 0.01, polycrystal states with the large particles forming grain boundaries for 
0.01 < c < 0.12, and glass states for larger c. For c > 0.1, the system is divided 
into small crystalline domains composed of the small particles and percolated amor- 
phous regions composed of the small and large particles. The areal fraction of the 
amorphous regions increases with increasing c. See Fig. 2 of Ref.l4 for the particle 
configurations for various c at cr2/o'i = 1-4, where a crystal state was realized for 
c = 0.02 with large particles being distributed as point defects. 

In this work, we carried out many simulation runs also for c > 0.5. For not very 
small c or 1 — c, the system is divided into crystalline regions composed of the large 
particles and amorphpus layers composed of the two species. The areal fraction of 
the layers is decreased for relatively small c or 1 — c. Figure 2 provides such an 
example for c = 0.9, where we display Dj{t) at t = 10^ in the left and the broken 
bonds and the displacement vectors 

Ar,{t,At) = rj{t)-rj{t- At), (16) 

in the time interval [10^, 4 x 10*^] with width Z\t = 3 x 10^ in the right. In writing the 
displacements we pick up only mobile particles with \Arj\ > 0.2. At this composition 
percolation of the amorphous layers is attained. However, for c = 0.95, the layers 
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Fig. 4. (color online) Disorder variables Dj{t) (upper plates), and broken bonds (x) and displace- 
ment vectors for mobile particles with \Arj\ > 0.2 (lower plates) for c = 0.3, 0.5, and 0.8 with 
7 = 0. For c = 0.3 and 0.5 the particles with Dj > 8 are written in black (upper plates). 



C=0.05 t=20000 X 30000 ❖ 40000 V 50000 o C=Q.3Q t-20000 X 30000 40000 V 50000 Q 




Fig. 5. (color online) Broken bonds and particle displacements with > 0.2 (arrows) in con- 

secutive four time intervals with width At = 10** for c = 0.05 (left) and 0.3 (right) with 7 = 
in the whole system. Shown is time evolution of the dynamic heterogeneity in polycrystal (left) 
and glass (right). 
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are broken into isolated clusters with various sizes, resulting in a crystal state with 
point defects (not shown here). In Fig. 2, the broken bonds are mostly created 
by the particle motions in the layers and are heterogeneously distributed. We also 
find coherent motions of crystalline domains surrounded by broken bonds at their 
grain boundaries. See the middle part at the bottom in the two panels for such an 
example. The patterns should be isotropic under the periodic boundary condition 
(on the average of many runs). Next, in Figs. 3 and 4, we show the same quantities 
for c = 0.05, 0.1, 0.15, 0.3, 0.5, and 0.8 to illustrate the crossover with increasing c. 
The system is in a well-defined polycrystal state for c = 0.05 and is highly disordered 
for c = 0.3 and 0.5. The structural heterogeneity in Dj and the dynamical one in 
the broken bonds and the displacements are both conspicuous for any c. Their 
correlation is obviously seen in the presence of small crystalline regions. This was 
illustrated in Ref. 14 in more expanded snapshots including the elusive case of glass. 

Furthermore, in Fig. 5, we show time evolution of the broken bond distributions 
in four consecutive time intervals with width 10"^ in polycrystal at c = 0.05 and in 
glass at c = 0.3. The structural relaxation time Tq, is of order 10^ at c = 0.05 and 10^ 
at c = 0. 3§ Here the clusters of the broken bonds mostly overlap or are adjacent to 
each other, resulting in aggregations of the broken bonds. Similar figures illustrating 
the heterogeneity evolution can be found in other papers 

§4. Numerical results at 7 = 10~^ 

We applied a simple shear flow with shear rate 7 = 10~^. To this end, we 
divided the system into three regions Hll'IHJ In the bulk region — 0.5L < x,y < 0.5L, 
we initially placed particles with = 9000. We added two boundary layers in the 
region — 0.6L < y < — 0.5L at the bottom and in the region 0.5L < y < 0.6L at the 
top. In each layer, N^, = 900 particles with the same composition and size ratio were 
initially placed. They were attached to the layer by the spring potential, 

Uj{r-R,) = ^K\r-R,\^, (17) 

where Rj = Rj{t) are pinning points in the boundary layers dependent on t. The 
spring constant was set equal to K = 20ea^^. These bound particles also interacted 
with the neighboring bound and unbound particles with the common Lennard-Jones 
potentials in Eq. ([s]) . The x component of Rj was moved as 

X,it)=X,{0)±hjt. (18) 

We imposed the periodic boundary condition in the x direction. When Xj{t) > L/2 
(or Xj{t) < —L/2) in the integration, Xj{t) was decreased (or increased) by L. The 
unbound particles rarely penetrated into the boundary layers deeper than ai. 

The unbound particles in the bulk regions obeyed the Newtonian equations of 
motion without thermostat. However, a Nose-Hoover thermostatic '1^ was attached 
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to each boundary layer independently. 
That is, the bound particles j G B were 
governed by 
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drj 



U -CBmj{rj-VB), (19) 
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Fig. 6. Averaged temperature T{y,t) vs y/L 
with 7 = 10~* for c = 0.05, which is higher 
in the middle than the boundary value 0.2 
only by a few %. 



where B denotes the top or the bottom. 
The boundary velocity i;g is (L7/2)e^ 
at the top and —{L^/2)ex at the bot- 
°2 0-4 06 torn, where e^, is the unit vector along 
the X axis. For 7 = 10~^ the boundary 
speed (~ 0.005) is much slower than the 
thermal velocity (~ 0.6). The two ther- 
mostat parameters, Cbot at the bottom 
and Ctop at the top, obeyed 



dt 



Cb 



1 



r, 



NH 
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(20) 



We set tnh = 0.304 and T = 0.2. For 7 = 10 the local temperature was nearly 
homogeneous in the bulk region during plastic flow, as demonstrated in Fig. 6 for 
c = 0.05. Setting At = 10^ and Ay = L/20, we define the averaged temperature by 



fiy,t) 



1 ft fy+Ay/2 r 

r J—r / dt' / dy' 
LAyAt Jt^^t Jy-Ay/2 J- 



y+Ay/2 ^ rL/2 

dy' / dx'eK{x',y',t'), 

Ay/2 J-L/2 



(21) 



where eK{x,y,t) = mj\Rj — '^yex\ 5{r — Rj)/2 is the kinetic energy density. 

We prepared the initial states as follows. First, we equilibrated the bulk and 
boundary regions independently without their mutual interactions. The particles in 
these regions interacted via the LJ potentials in Eq.Q in a liquid state at T = 2 in 
a time interval of 10^ under the periodic boundary condition along the x and y axes. 
Second, we quenched the system to T = 0.2 and further equilibrated it for a time 
interval of 10'^. Afterwards, we chose the particle positions in the boundary layers 
as the initial pinning points Rj{0) and introduced the spring potential of the bound 
particles and the LJ potentials between the bound and unbound particles. Third, 
we further waited for time interval of 5 x 10^ until we detected no appreciable time 
evolution in various quantities. After this second low-temperature equilibration, we 
applied a shear flow with rate 7 = 10~^ by sliding the pinning positions in the 
boundary layers with velocities ±7^/2. We set t = at the application of shear. 

In Fig. 7, we show the broken bonds in five consecutive time intervals with width 
5000 = 0.5/7 ™ polycrystal at c = 0.05. Lines of the broken bonds represent plastic 
events created by dislocation glidin^lD or grain boundary slidinglHf The former 
deformation mode is dominant at this small c, while the latter becomes increasing 
important with increasing cH^'llll' We can see that most of the slip lines are nearly 
parallel to the x or y axis. This is because the elastic energy of a slip is minimum in 
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Fig. 7. Broken bonds in five consecutive time intervals with widtii At = 5000 = O.5/7 in shear 
7 — 10~* in polycrystal with c = 0.05. They form slip lines connecting grains. The top and 
bottom boundary layers are without broken bonds. Right bottom: averaged velocity Vx{y,t) 
defined by Eqs. (12) and (13) from the same data for these time intervals, which greatly changes 
in time and much deviates from the linear profile. 



these directions in simple shear strain. ' In the first three time intervals a thick band 
region extends over the system along the x axis in the upper part, but in the last 
two time intervals such large-scale bands are nonexistent. The right bottom panel 
displays the average velocity Vx{y, t) defined in Eqs. (12) and (13) with At = 5x 10^, 
where the first three curves largely deviate from the linear profile jy but the last 
two curves are nearly linear. 

In Fig. 8, we show the same quantities in five consecutive time intervals with 
width 2500 = O.25/7 in glass at c = 0.3. Also in this case, the plastic deformations 
are highly heterogeneous, which still tend to be nearly parallel to the x or y axis. The 
chains of the broken bonds are shorter than in Fig. 7. That is, the slip lines do not 
much exceed the molecular size and the broken bonds are more randomly scattered 
than in polycrystal. As a result, the deviation of the averaged velocity Vx{y, t) from 
the linear form becomes small for At > I/7, while it is large for At < O.I/7 (see Fig. 
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Fig. 8. Broken bonds in five consecutive time intervals with width At — 2500 — 0.25/7 in shear 
7 — 10~* in glass with c = 0.3. Slip lines are short. Right bottom: averaged velocity Vx{y,t) 
from the same data for these time intervals. 



9). In their long time simulation of sheared 2D glass, Furukawa et al^i realized 
transient shear bands equally in the flow and velocity-gradient directions with the 
Lees-Edwards boundary condition, while those parallel to the x axis more easily 
extend upto the system length in our simulation with the boundary layers. 

Figure 9 illustrates the time-evolution of the dynamic heterogeneity in three 
consecutive time intervals with width 400 = O.O4/7 in polycrystal at c = 0.05 and 
in glass at c = 0.3. We display the broken bonds in the left panels and the averaged 
velocity Vx{y,t) in the right panels. The width At here is made much shorter than 
in Figs. 7 and 8 to analyze the early stage of plastic deformations. Note that the 
transverse sound velocity c± is of order 5 (in units of cr/r) and the acoustic traversal 
time L/c± is of order 20 in our system. (Plastic deformation dynamics on shoter time 
scales has been analyzed in our previous work.^-^*) In the upper panels at c = 0.05, 
large-scale plastic deformations multiply occur in the same area, which extends over 
the system to form a shear band. On the other hand, in the lower panels at c = 0.3, 
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Fig. 9. (color online) Left: time evolution of bond breakage in three consecutive time intervals with 
width At = 400 = 0.04/7 in shear 7 = 10"". Right: averaged velocity Vx{y,t) for these time 
intervals much deviating from the linear profile. The system is in polycrystal with c = 0.05 (top) 
and in glass with c = 0.2 (bottom). Repeated occurrences of plastic events result in transient 
shear bands. 



a large shear band is seen in the lower part in the first time interval, but large plastic 
deformations occur in the vertical direction in the second time interval and in the 
upper part in the horizontal direction in the third time interval. The life time of the 
shear bands is shorter in glass than in polycrystal.^^' In addition, the broken bond 
numbers ANf, in the whole system are 649, 551, and 614 for c = 0.05 and are 870, 
589, and 761 for c = 0.3 in these time intervals in the chronological order. The time 
average of ANb is of order NjAt{= 360 here), which should be the case for any c as 
discussed at the end of Subsec.2.3llll' 
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§5. Summary and remarks 

We have examined the jammed particle configurations and the dynamic het- 
erogeneity in 2D. Visualization of the disorder variable Dj gives information of the 
structural heterogeneity, while that of the bond breakage and the particle displace- 
ments discloses the presence of the dynamic heterogeneity. We have varied the 
composition c for cr2/cJi = 1.4 and T = 0.2e/kB- Summarizing our main results, we 
give some remarks below. 

The particle configurations are unique in polycrystal and glass which are realized 
for not very small c or 1 — c in our 2D model system. As in Figs. 2-4, small crystalline 
regions composed of one species of the particles are enclosed by percolated amorphous 
regions composed of the two species, where the amorphous regions form layers for 
relatively small c or 1 — c. For very small c or 1 — c (not shown in this work), the 
layers break into small pieces to form a crystal with point defects. In rheology, these 
amorphous layers can serve as a lubricant in plastic flowjlll' within which the bond 
breakage occurs, reducing the viscosity than in crystal states. 

In glass without shear, the particle configuration changes are thermally induced 
in the form of chains of broken bon or stringlike displacements;^^' They 

accumulate to form the dynamic heterogeneity on long times characterized by the 
correlation length ^ in Eq.(l)ll2)'[i§'ll3 In Figs. 2-4, we have added more evidences 
for these dynamical processes. To show the relationship between the structural and 
dynamical heterogeneities, we have presented the snapshots of the disorder variable 
Dj and those of the bond breakage and the displacements (see Fig. 2 in Ref.l4 
also). In Fig. 5, we have visualized the time evolution of the dynamic heterogeneity 
on the time scale of 10^, which demonstrates correlated occurrences of mesoscopic 
structural changes. This tendency is related to the mesoscopic heterogeneity in the 
elastic modulilH2l'll3 In addition, in Fig. 4, the structural relaxation is even slower 
in polycrystal than in glass, where the particles around the grain boundaries (in the 
amorphous layers) are relatively mobile than those within the crystalline grainsEH' 

In shear flow with 7 = 10~^, we have visualized large-scale heterogeneity in the 
bond breakage in Figs. 6-8. The averaged velocity Vx{y,t) deviates from the linear 
profile with formation of shear bands. Here At = 0.5/ j for c = 0.05 in Fig. 7, 
At = 0.25/7 for c = 0.3 in Fig. 8, and At = O.O4/7 for c = 0.05 and 0.3 in Fig. 9. 
These figures illustrate how the plastic deformations under shear evolve on various 
time scales. For any c, plastic deformations extend over longer distances nearly in 
the X or y axis as in Fig. 9. Such plastic events often accumulate to form transient 
shear bands as in Figs. 7-8. They occur to release the elastic energy at high strain 
as discussed in Sec. I. See Refs. 27 and 28 for more discussions on the anisotropic, 
hierarchical dynamics under shear. In the early work by Yamamoto and one of the 
present authors jHI the shear- induced structural change was analyzed on the basis 
of Eq.(2) and the broken-bond structure factor Sb{k) in Eq.(l) was calculated for 
At = O.OSrft = 0.5ro without and with shear, but the shear band formation on large 
scales was beyond its scope. 

A variety of complex problems of binary particle systems remain mostly un- 
explored. We mention some of them, (i) When the size ratio a^jox is increased 
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from unity, intriguing crossovers are expected in the particle configurations and the 
dynamics. For example, at c = 0.5, proliferation of defects occurs abruptly around 
a^jox — 1.20''[llf (ii) Furthermore, with changing the pair potentials, tendency 
to phase sparation can be enhanced, where small crystalline regions should become 
more disinct or nucleation of crystal domains should become realizable in an amor- 
phous matrix, (iii) Some essential aspects of the glass problem should be common 
in 2D and 3D. However, in 3D, the problem is much more complicated and visual- 
ization of configuration changes is more difficult. We have not yet understood the 
differences in 2D and 3D in depth, (iv) In addition, the melting and crystallization 
in binary mixtures have not yet been well investigated, where essentially different 
pictures arise in 2D and 3D. We have recently recognized surprising complexity of 
the 2D melting even in one component systems!^ That is, marked heterogeneities 
emerge both in the structural disorder and in the dynamics in the hexatic phase jlSl 
where mesoscopic liquidlike and crystalline regions coexist as thermal fluctuations 
without distinct interfaces. Understanding the dynamics of 2D melting is still at the 
beginning despite numerous papers on static properties. 
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